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Abstract 



We use the BFKL equation to calculate the rate of deep inelastic scattering events 
containing two forward jets (adjacent to the proton remnants) at HERA. We compare the 
production of two forward jets with that of only one forward jet (the "Mueller" process). 
We obtain a stable prediction for this two to one jet ratio, which may serve as a measure 
of the BFKL vertex function. 



I. Introduction 

One of the most striking discoveries at the electron-proton coUider, HERA, at DESY, is the 
steep rise of the proton structure function F2{x,Q'^) with decreasing Bjorken x = Q'^/2p ■ q, 
where = — g^. Here p is the four-momentum of the proton and q is the four-momentum 
transfer of the electron probe. The behaviour of the structure function at small x is driven by the 
gluon through the process g qq. The behaviour of the gluon distribution at small x is itself 
predicted from perturbative QCD via the Balitskij-Fadin-Kuraev-Lipatov (BFKL) equation, 
which resums the leading [as ln(l/x)]" contributions. This predicts a characteristic x"'^ singular 
behaviour in the small x regime which results from the summation of soft multi-gluon emissions, 
which effectively correspond to the sum of gluon ladders (Fig. 1) where the transverse momenta 
kn of the emitted gluons are unordered along the chain. For fixed as the BFKL exponent A = 
as4 In 2, where as = Sas/vr. It is this increase in the gluon distribution with decreasing x which 
produces the corresponding rise of the structure function. Of course there are uncertainties due 
to subleading corrections, which will reduce the value of the exponent A, and from the treatment 
of the infrared region in the numerical solution of the BFKL equation. A recent description of 
F2 along these lines, but also including the effects of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution [Q, can be found in ref. 

However, the situation is more complicated. DGLAP evolution, which resums the as In 
contribution can on its own produce an increase in F2 with decreasing x, and with a suitable 
choice of input distributions at scale Ql can give an excellent description of F2 at small x. 
This makes it difficult to distinguish whether the steep rise observed for the structure function 
has a major component from ln(l/x) effects, or is simply due to InQ^ evolution from a non- 
perturbative input. The observable F2 is just too inclusive to act as a discriminator of the 
underlying small x dynamics. 

In 1990 Mueller p proposed that the less inclusive quantity of DIS plus an identified 
energetic forward jet may serve as a useful probe of small x dynamics. The idea is to isolate 
DIS (x, Q"^) events containing a jet {xj, k^j) with x <^ Xj and k^j ~ Q^, see Fig. 1. Here kxj and 
Xjp are the transverse and longitudinal momentum components of the jet. We can thus expect 
a {x/xj)~^ behaviour in a kinematic region free of DGLAP evolution^. Moreover, for large 
Xj ~ 0.1, the parton distributions of the proton are well known and so the Mueller proposal 
removes the uncertainty arising from the introduction of a non-perturbative input. That is we 
study DIS off a known parton, rather than off the proton. Recently data for the process have 
become available and the comparison with BFKL predictions is encouraging 

The observed DIS + forward jet data sample contains a fraction of events in which two jets 
are identified, each with pt greater than the experimental cut. Our concern here is to extend 
the BFKL formalism to study and to estimate the rate for this process. We also present the 
predictions as the ratio of 1 to 2 forward jet production. 

^The [as In Q^]" contributions which are resummed by DGLAP come from the region in which the transverse 
momenta of the emitted gluons are strongly ordered along the chain, that is 3> 

■> ...klj. DGLAP 

evolution is thus neutralized by the choice ~ k'^y 
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The structure of the paper is as follows. In Sec. II we use the data for DIS + one forward 
jet to determine the normalisation of our numerical solution of the BFKL equation. In Sec. Ill 
we discuss the extension of the BFKL formalism to include an extra identified gluon jet, and 
then in Sec. IV we present results for the cross section of DIS + two identified forward jets. 
Our conclusions are given in Sec. V. 



II. DIS + one forward jet 

The differential cross section for deep inelastic scattering (x, Q^) containing a forward jet 
{xj, k^j) is of the form 



da 



Ana' 



dxdQ'^dxjdkj.j xQ'^ 



+ w 



dxj dk'^j ^ dxjdkj 



(1) 



where y = Q'^/sx with y/s equal to the centre-of-mass energy of the electron-proton collision. 
The differential structure functions 



OF 



3as{k^ 



doc j dk^j 



irkj^j 



''^^jfa{xj,krpj)\ ^i(x/xj, kjnpQ ) 



(2) 



where the $j represent the photon-gluon process shown by the blob in Fig. 2(a). The factor k^j 
arises from the gluon propagators. We are interested in events with small x/xj so Xj should be 
taken as large as is experimentally feasible. Eq. (2) assumes strong ordering of the longitudinal 
momentum at the gluon-parton a vertex, so that Xj of the exchanged parton, which occurs in 
the parton density fa{xj, k^j), is to a good approximation that of the outgoing jet. Assuming 
t channel pole dominance we have 



E /a = 9 + liq + q) 



(3) 



where g, q and q are the gluon, quark and antiquark densities respectively. To include the 
running of as it is convenient to work in terms of the functions [H 



H,{z,e,) = as{kl)^^{z,kl,Q'') 
with i = 2,T. These functions satisfy the BFKL equations 



(4) 



(5) 



1 dz' f d^T 



+ -.(4) / ^ / 



rrq^ 



'^H,{z',k'S,Q^) - H,{z',kl,Q^)Q{kl-ql) 



L.'2 

ri/rp 



where k'^ = [kx + Note that the g|. — > singularity is removed by the cancellation 

between the real and virtual contributions. Due to the strong ordering at the gaa vertex the 
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transverse momentum of the BFKL gluon is that of the forward jet, that is ~ k^j. In 
eq.(5) we impose a cut-off k'^ {k"^) > k^ (with k^ ~ IGeV^) in order to avoid the Landau pole 
singularity of as{k'^). The driving terms = as{kj)^f^ are given by the quark box (plus 
the crossed box) contributions to the 7(7 fusion process 



K K — fcr\^ o/l 1 



where k' = k — {1 — P)kT and 

D2 = (k - krf + p{l - P)Q^ + ml 
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(7) 



For small z, eq.(5) generates the characteristic z^'^ behaviour with A = 41n2 for fixed 
coupling as- DIS + forward jet events have been discussed within the BFKL formalism with 
fixed coupling 0. 

Our objective is to calculate the ratio of DIS events containing two, as compared to one, 
forward jets. Since the contributions from the input largely cancel, for simplicity we do 
not include charm in the equations for the quark box and we assume massless quarks. In the 
z' integration in (H) we impose a cut z < Zq 1, such that Xj > x. This determines the 
point, Zq, at which the evolution in ln(xj/x) starts, that is the point at which the gluon chain 
is matched onto the quark box . Essentially, we take Zq as a free parameter which is chosen 
to normalise the prediction for the DIS + one forward jet cross section (as a function of x) to 
the measured values. We find zo ~ 0.5. If we were to enhance the driving terms by 20%, say, 
to allow for a charm contribution, then the starting input value zq would be about 0.1. 

To make the comparison between the prediction and the data we must impose the exper- 
imental kinematic cuts used in the measurement of the DIS + jet events. We adopt the cuts 
used by the HI collaboration 0]. That is the forward jet is required to lie in the region 

7° < Oj < 20°, Ej > 28.7GeV, krj > 3.5GeV (8) 

whereas the outgoing electron is constrained to the domain 

160° < < 173°, > llGeV, y > 0.1 (9) 
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in the HERA frame. Finally HI require that 



|g2 < < (10) 

so as to minimize the DGLAP InQ^ effects. 

We use these cuts to calculate the DIS + forward jet cross section and compare with the HI 
data in Fig. 3. We find that the choice zq = 0.5 leads to a good description of the data. It is 
relevant for our study of the 2 forward jet rate to note the effect of the upper phase space cut 
on the variable y. The solid and dashed histograms in Fig. 3 correspond to using y < 0.5 and 
y < I respectively. We see that the predictions at small x are sensitive to this cut. In fact if 
the next smallest x bin, Ax = (0.0001, 0.0005) had been shown we would have seen a turnover 
due to the depletion of events caused by the lack of phase space. 



III. DIS + two forward jets 

Here we extend the BFKL formalism so as to be able to predict the rate for the production 
of an extra jet in the forward direction. The process 

7*+P^Ji+j2 + X (11) 

is shown in Fig. 2(b). We require that both jets are resolved and have transverse momentum 
greater than some minimum cut, k^j^ > ^"^ with i = 1,2. The differential cross section for this 
process is given by 

da 



dxdQ"^ dxji dk^j^ dxj2 dk'^j2 



(12) 



47ra 
1^ 
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dxji dkj^j^ dxj2 dkq.j2 ^ 9xji dkj,j^ dxj2 dkj,j2 



where the differential structure functions 

-Xj2 



^ dxji dkj.j-^ dxj2 dkj'j2 

(13) 



where the variables are shown in Fig. 2(b). The variable cj) is the azimuthal angle between the 
two forward jets. Because 

kl = {kTji + kTj2? 

the integral over is non-trivial. The general structure of formula (p!3D for the production 
of two forward jets is readily understood as an extension of the single jet formula, (0). The 



BFKL functions $i describe the gluon radiation in the upper part of the gluon emission chain, 
which is shown as a blob in the Fig. 2(b). These functions are obtained by solving the BFKL 
equation exactly as described in Sec. II, and are indeed normalised so as to describe the DIS + 
one forward jet data. Since both of the resolved jets are required to be in the forward region 
Xj2 ^ Xji ~ 0(1), we do not enter the strongly ordered configuration Xj2 <C Xji ~ 0(1) and 
so we can neglect the effect of soft (BFKL) gluon radiation emitted in the rapidity interval 
between the two jets. The remaining noteworthy feature of (|13]) is the presence of the BFKL 
vertex function, k'^ji/k^j2 controlling the emission of jet 2. 

When the second jet is produced in the central region (that is when Xj2 <^ xji) then the 
gluon radiation in the rapidity interval between the two jets can no longer be neglected. We 
would need to use the formalism developed in ref. [^. 

IV. Predictions of the DIS + two forward jet rate 

We calculate the cross section for DIS + two forward jets from eqs.([I^) and (0), using the 
kinematic cuts listed in eqs.(|8|-|T0D. We present the cross section in Fig. 4. To gain insight into 
the effect of the cuts we also include the results obtained using the enlarged domain 

|Q2 < A;2^. < 5Q2 (^4) 

for each of the two jets. For the experimental identification of two jets we require a minimum 
angular separation between the jets. We therefore impose a cut on the [r], 0) space of the two 
partons producing the forward jets. For illustration we show the results for 

^/iA^i)^TJA^ > R (15) 

for four different values of R, namely R = 0,1, 1.7 and 2, where Arj and are the differences 
in the pseudorapidities and azimuthal angles of the two parton jets. 

From Fig. 4 we see that the two jet cross section has a turnover for the smallest x bin for the 
lowest set of curves , which correspond to the cut k^j/Q'^ < 2. The effect is not apparent when 
we relax the constraint to k^^/Q'^ < 5. Also we note that the cross section is appreciably larger 
for the higher kinematic cut, k^j/Q'^ < 5. First we quantify the increase. When we change the 
krj/Q"^ cut from 2 to 5 the cross section for DIS + single forward jet increases by a factor of 
about 2.5 for the smallest x bin that is shown, Aa; = (0.0005, 0.001), but increases only by a 
factor of about 1.5 for the largest x bin. Ax = (0.003, 0.0035). Similarly, the cross section for 
DIS + two forward jets increases by factors of 5.5 and 2.5 for these two Ax bins respectively. 
The explanation for the large increase of cross section as we relax the upper limit on k^j from 
2Q^ to 5Q^ can be found by inspecting Fig. 5. 

Fig. 5 shows the differential cross section for DIS + two forward jets, for the x bin Ax = 
(0.001, 0.0015), as a function of the transverse momentum of the jets, which are taken to have 
equal magnitude k^j-^ = k^j2- We take jet resolution cut, k'^j > /i^, to be /i = 3.5GeV. The 
curves on Fig. 5 correspond to the differential cross section for different values of Q^. When 
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-0.001 
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39.6 


0.001 


- 0.0015 


6.0 
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19.4 


25.2 


0.0015 


-0.002 


4.4 


8.6 


12.3 


15.3 


0.002 


- 0.0025 


3.2 


5.9 


8.0 


9.7 


0.0025 


- 0.003 


2.3 


4.1 


5.4 


6.4 


0.003 


- 0.0035 


1.7 


2.8 


3.7 


4.3 



Table 1: The dependence of the DIS + two forward jet cross section (in pb) on the kinematic 
constraint ^Q^ < k'^j < nQ"^ for various values of n. We show the dependence for different 
intervals of x. We take i? = 0, y <1 and ji = 3.5GeV. 

we make the cut k^j/Q'^ < 2, we remove a substantial contribution to the cross section from 
the low region. The effect is enhanced by the jet resolution cut k^j > fi^, which requires 
the cross section to be zero when fi'^/Q'^ < 2, as indicated by the vertical arrows on Fig. 5, 
which occur above the lower bound k'^^/Q'^ = | for both = 15 and 12GeV^ if we take 
H = 3.5GeV. Thus the resolution constraint has significant impact at the lower values of Q^. 
The X dependence of the cross section for different kinematic constraints is shown in Table 1. 

The turnover that is apparent at small x in the lower set of histograms in Fig. 4 is simply 
due to the lack of phase space. The reason the cross section predictions for the smaller x bins 
are particularly sensitive to phase space restrictions, is that such cuts kill the important low 

contributions, see Fig. 5. 

The BFKL functions are common to the calculation of both the single and two forward 



jet rates, see eqs.(|D and (|T^). Thus the ambiguities in the calculation should be reduced in the 
prediction of the ratio cr2}et/o'i}ct- The predictions are shown in Fig. 6, and also Fig. 7 which 
shows the dependence of the ratio on the jet resolution parameter /x (but using larger Ax bins). 
As anticipated from the previous discussion we see that the ratio increases significantly when 
we relax the cut from k^j/Q'^ < 2 to k'^j/Q'^ < 5. We see that the ratio shows an increase as 
we move from lower to higher values of x. The reason is that the ordering Xj2 < Xji of the two 
emitted jets means that the reduction of the cross section due to decreasing evolution length 
is less for 2 forward jets than for single jet production. The ratio is insensitive to the choice of 
the upper limit for y. This is illustrated in Table 2 for two different jet separation cuts. 

We conclude that the ratio is an observable which can give information on the underlying 
small X dynamics, and in particular can act as a measure of the BFKL vertex function which 
occurs in (|T3p. The first experimental estimates of the ratio have been made [^], but so far 
only to give an indication of the size of the effect. A value of about 4% is quoted for the ratio 
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Table 2: The dependence of the 2:1 forward jet ratio on the kinematical [k^j/Q"^) and phase 
space (y) cuts. The cross section is shown for the x bin Aa; = 0.0005 — 0.002. fi is the jet 
resolution parameter, k^j > /x^, and R the jet separation parameter of eg. ([T5|) . 



of 2-jet/l-jet events for a slightly different choice of kinematic variables than those of (^)-([T^) 
above, namely for the domain 

6° < 9j < 20°, Xj > 0.025, krj > 5GeV, 

160° <9e< 173°, Ee > 12GeV, 0.1 < y < 1, 

0.0002 <x< 0.002, < k'^T < 4<5^- 

2 ■' 

For this kinematic region we predict a 2-jet/l-jet ratio of 6.6% at the parton level. 

A characteristic feature of the effects of the BFKL chain is the angular decorrelation of the 



jets [10|. Due to the lack of /cT-ordering the ku gluon of Fig. 2(b) can bring significant transverse 
momentum into the two jet system and hence considerably broaden the "back-to-back" peak 
in the azimuthal distribution A0 = 0i — 02- Sample distributions are shown in Fig. 8. 

In order to remove the infrared k^ ^ infinities in the BFKL equation we impose the cut 
k^ = {krji + kTj2f > IGeV^. This, together with the ambiguity due to hadronization, means 
that the distribution near the back-to-back configuration cannot be predicted. Rather it is the 
tails of the distribution which will characterize the lack of kx ordering. Unfortunately the jet 
separation cut, which experimentally is likely to be i? ~ 1.7, will effectively remove the whole 
of the tail of the distribution, see Fig. 8. The azimuthal decorrelation is thus unlikely to be a 
way of identifying the underlying small x dynamics. 



V. Conclusions 

The measurement of DIS scattering events (x, Q^) containing a very forward energetic jet 
{xj, k^j), with Xj ^ X and k^j ~ Q^, has long been advocated as a favourable way of investi- 
gating the dynamics which underlie small x physics. The forward region is defined by cuts of 
the type given in eq.(H) and (H), and jet events are collected which typically lie in the region 
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subject also to the jet resolution cut, k'^j > /i^. To obtain sufficient statistics the ffist measure- 
ments take /X = 3.5GeV, but ideally as the integrated luminosity at HERA improves, it will be 
preferable to obtain a DIS data sample restricted to jets with higher transverse momentum. It 
is already known that the DIS + forward jet data sample contains a small fraction of events 
with two identified jets The main purpose of this paper is to predict the fraction, and to 
study the properties, of these events. Our results are shown by the lower set of histograms in 
Fig. 6 for yU = 3.5GeV, and in Fig. 7(a) for three different values of fi. We see that the fraction 
of 2 jet events for R = 1.7 is about 2.5% for the x interval (0.5 — 1) x 10~'^ rising to about 5% 
for X = (2.5 — 3) X 10^^, see Fig. 6. The fraction is rather insensitive (i) to the value of the 
jet resolution parameter /i, see Fig. 7(a), and (ii) to ambiguities in the function $j describing 
the BFKL gluon chain since it is common to both the 1 jet and 2 jet predictions and tends to 
cancel in the ratio. The experimental conffimation of the predicted 2 jet fractions will therefore 
serve as a check on the BFKL vertex function which occurs in (p!3|). 

The domain ([l6| ) is chosen such that k^j ~ so as to suppress DGLAP gluon emission. 
We also presented results for the larger domain < k^j < 5Q^ simply to gain insight into 
the behaviour of the DIS + two forward jet cross section. 
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Figure Captions 



Fig. 1 Diagrammatic representation of a deep inelastic event containing an identified forward 
jet with longitudinal and transverse momentum components Xjp and kTj respectively. 
The photon scatters from the gluon chain (via the quark 'box') described by the BFKL 
equation which resums the (soft) gluon emissions. Parton a can be either a gluon or a 
quark. 

Fig. 2 The diagrams giving the cross section for (a) DIS + one forward jet and (b) DIS + two 
forward jets. In each case the upper blob represents the BFKL functions $j (which in the 
absence of the BFKL gluon emissions would, to lowest order, be simply given by the quark 
box and crossed box, ^f'^)- The dot represents the BFKL vertex function "measured" by 
the 2 jet/1 jet ratio. 

Fig. 3 The DIS + forward jet cross section compared with recent HI data @]. The phase space 
cuts are given in eqs.(|8|- |T0|) . The resummed BFKL kernel is normalized by the choice 
zo = 0.5. The solid (dashed) histograms correspond to y < 0.5 and y < I respectively. 

Fig. 4 The DIS + two forward jet cross section (72 in the same x bins as for the DIS + forward 
jet results. We take y < 0.5. a2 is plotted for krj > 3.5GeV for two different choices of 
the upper krj cut: k^j/Q"^ < 2 and 5. The results are shown for four different cuts on 
the minimum separation between the two forward jets: R = 0,1, 1.7 and 2. The smaller 
the separation cut, the higher the cross section. 

Fig. 5 The differential cross section for DIS + two forward jets as a function of k'^j/Q'^, for 
different values of Q^. We take k'^j^ = k'^-2, k\,-IQ'^ > |, i? = and the jet resolution 
cut off /i = 3.5GeV. The cross section ([T^) is integrated over Xj\ and Xji and is shown 
for the X bin Ax = (0.001, 0.0015). The lower cut-offs are indicated by vertical arrows 
for the different values of . 

Fig. 6 The same as Fig. 4, but showing the ratio of the two/one forward jet cross section. 

Fig. 7 The two/one forward jet ratio predicted for DIS events shown for three different values 
of the jet resolution parameter /i (that is > jjp'), and for two upper limits of k^j/Q"^. 
The two jets are taken to have the same k^j. The five different histograms correspond, 
in descending order, to i? = 0, 1, 1.5, 1.7 and 2, where the jet separation parameter R 



is defined in (|T5|). 



Fig. 8 The cross section as a function of the azimuthal separation, A0 = 01 — 02, of the two 
forward jets. The distribution is shown for two choices of k^^/Q"^ and for two different x 
bins. We take y < 0.5, and we show the distributions for i? = 0, 1 and 1.7, where the jet 
separation parameter R is defined in (Kf). 
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